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We present an operational approach to address the in-medium behavior of charmonium and ana- 
lyze the reliability of maximum entropy method (MEM). We study the dependences of the ratio 
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tinuum's threshold. Furthermore, we discuss the issue of the default model dependence of the 
spectral function obtained from MEM. 
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1. Introduction 

Matsui and Satz [|l|] proposed long time ago the melting of charmonium states due to color 
screening can signal the formation of Quark Gluon Plasma in heavy ion collisions. The suppression 
of the electrons from J/y is observed both at RHIC and SPS [0], however, its interpretation is not 
quite understood. This phenomenon has been studied quite extensively in potential models y], 
based either on models or finite temperature lattice QCD results for the heavy quark potential in a 
non-relativistic Schrodinger equation. Depending on the form of the heavy quark potential used, 
the dissociation temperature of charmonium can be ranged from shortly above T c up to values 
similar to those obtained in lattice QCD [Q]. With the lattice QCD approach, the properties of 
the charmonium, which can be directly seen from the spectral function, is enclosed in the lattice 
QCD calculated euclidean time correlation functions. To extract the spectral functions from the 
correlation function, Maximum Entropy Method ^ is normally used. 

Here we contribute an operational approach to address the in-medium behavior of charmonium 
and address the issue of the default model dependence of the spectral function obtained from MEM. 



2. Lattice correlators and spectral functions 

We look into the momentum-projected Matsubara correlators 

G(T,p,r)=^e^<7 fl (T,x)4(0,0) > T , (2.1) 

x 

where Jh is a suitable mesonic operator, p is the spatial momentum, T is the temperature of the 
gluonic plasma and the Euclidean time z € [0,1/T). Through analytic calculation, the Matsubara 
correlator can be related to the hadronic spectral function as the following: 



G H (z,p)= d<oo H (<o,p,T) \ (2-2) 



COSfl((0(X— 

sinh(- 2T , 

Extracting the spectral function at finite temperature lattice QCD is hampered mainly by two 
issues: the physical extent of time is restricted by the temperature, X < l/T, and the finite number 
of correlator points 1 making the inversion of Eq. 2.2 ill-posed. 



3. Charmonium correlators 



First, we analyze the sensitivity of the correlators to the spectral function by using the two 
following references correlators: 

poo poo 

G (t,T)=J o dcoa(co,T = 0)K(T,T), G free (z,T) = J dcoa free ((0,T)K(x,T), (3.1) 

where Gq(x,T) is so called "reconstructed" correlator, and Gf ree (z,T) is the free correlator at 
finite T. Gq(x,T) and Gf ree (z, T) show the behavior if the spectral function at temperature T were 

'what's more, the correlator points are not precise but with statistical errors 
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Figure 1: Ratio of correlator G(t, T) to Gq(t, T) versus tT: with different width of the resonance (left) and 
different threshold of the continuum (right). 

identical to that at T=0 and at free limit, respectively. The ratios of the finite T correlators to these 
two references are: 

G(r,T) G(r,T) 
G {Z,T) G free {T,T) 
then /?o ~ 1 and R f ree « 1 indicate the bound state still exist and already melt, respectively. 

We consider the spectral function as a combination of the resonance and the continuum: a = 
o r + a c . For the spectral function of the resonance, the following form is taken for T=0: G r ((0, T = 
0) = 8((0 — M), where M donates the mass of resonance. And at finite temperature T, we take 
the spectral function of the resonance to have the relativistic Breit-Wigner form of o r {(0,y) = 
}, where / is the width of the resonance, and N(y) is the normalization factor 
to maintain the relativistic Breit-Wigner the same strength as the delta function. For the continuum 
part of the spectral function, we take the formula of o c = ■^2(D 2 tanh(-^)yJ\ — (^) 2 (2 + (^) 2 ), 
where s is threshold of the continuum, for T=0, s = so = 4.5 GeV, for the free case, s=2m (m is 
mass of the quark), for finite T, s is T dependent. 

The ratio Ro is shown in the Fig. 1, in which the left plot is with different values of the 
resonance's width and the right plot is with different values of the continuum's threshold. We can 
see that both increasing the resonance's width and decreasing the continuum's threshold can make 
Ro go farther away from the unity. The influence of the resonance, with the width of 0.9 GeV only 
making a deviation of 9% at the symmetry point, is much smaller than that of the continuum, with 
threshold being 0.8 GeV smaller than so=4.5 GeV making a difference of 20%. The ratio Rf ree with 
different values of the width of resonance and different values of the threshold of the continuum is 
shown in the left and right plot of Fig. 2, respectively. Similar to Rq, the influence of the resonance 
is smaller than that of the continuum. 



4. Reliability of MEM 

After checking the charmonium on the correlator level, we're going to the spectral function 
level. The normal technique to extract the spectral functions from the correlator is the MEM, 
by maximizing a function Q(a;a) = a5[cr] — L[a]. L[a] is the usual likelihood function and 
minimized in the standard % 2 fit. The Shannon- Jayes entropy S[a] is defined as 

r°° a(co) 
S[a}= (tocy 0) )-m {a )-a {(o)log , (4.1) 

Jo mica) 
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Figure 2: Ratio of correlator G(t, T) to Gf ree (x, T) versus tT: with different width of the resonance (left) 
and different threshold of the continuum (right). 

where m(to) is the default model (DM hereafter) and should be given as a plausible form of a(co). 
a is a real and positive parameter which controls the relative weight of the entropy S and the 
likelihood function L. The final output o out is determined from a weighted average over a, G out ~ 
/ Aao a {(o)P[a\Dm], where the most probable spectral function o a {(Q) for given a is obtained by 
maximizing the Q and P[a\Dm] is the weight factor. The DM is very important as it strongly 




Figure 3: The DM dependence of spectral function obtained from MEM. 

affects the output of the MEM when the quality of the data is not sufficient. As we're focusing 
on the modification of the ground state of the spectral function, the nature choice of the DM is the 
asymptotic behavior of the spectral function at large CO in the free limit. Fig. 3 shows the outputs 
of the MEM when using different default models (m(co) = F * ^i® 2 )- The black line is the input 
spectral function (one resonance plus continuum) and the mock data is obtained by adding random 
Gaussian noises. All the three default models reproduce the location of the resonance well, in 
which the one with F=l (has the same large CO behavior as the input spectral function) gives the 
most reliable image. For the one with F=3, the output gets wiggled after the resonance, which is 
normally considered as "lattice artifacts" but could also be the "MEM artifacts" At this point, 
it could be better to use the free lattice spectral function [^] as the DM in practice. 

However, due to the lattice cutoff, such an asymptotic behavior (like DM with F=l) is not so 
obvious and can not be obtained directly in the lattice simulation. It could be helpful to look into 
the weight factor distributions and the correlators calculated from the default models, which are 
shown in Fig. 4. When one puts some physical prior information into the DM, it could be better 
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Figure 4: Weight factor distribution (left, for the visibility, the green line and the blue line are vertically 
multiplied by 10 and 10 3 , respectively), correlators calculated from the DM and the input correlators (right). 

that the correlators calculated from the DM is somehow comparable to the lattice correlator data 
(and consequently the larger peak location or amplitude of the weight factor function, see left panel 
of Fig. 4) rather than some orders of differences (see right panel of Fig. 4). And one also has to 
keep in mind that, no matter what kind of DM used, the correlators, that calculated from the output 
spectral functions that obtained from the MEM, can always reproduce the lattice correlator data 
within the errors. This essentially accents the importance of the prior knowledge of the spectral 
function to put into the DM and a careful analysis of the DM dependence. 

5. Summary 

Within current scenario of the spectral function, the correlator is more sensitive to the change 
of the continuum part than the resonance part, which makes the exploration of the resonance's 
properties more difficult. For the frequently used method MEM, it may also produce some arti- 
facts. As MEM can always reproduce the lattice correlator data, it's very important to put as much 
physical information as possible into the DM. When one is suspicious about some parts of the 
spectral function from the MEM, it could be helpful to look into the correlators calculated from the 
default models and the weight factor distributions. In this analysis the zero mode contribution [[7|] 
is not included, which requires further research. 
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